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ABSTRACT 

Context. OB stars are important constituents for the ecology of the Universe, and there are only a few studies on their pulsational 
properties detailed enough to provide important feedback on current evolutionary models. 

Aims. Our goal is to analyse and interpret the behaviour present in the CoRoT light curve of the B0.5IV star HD51756 observed 
during the second long run of the space mission, and to determine the fundamental stellar parameters from ground-based spectroscopy 
gathered with the Coralie and Harps instruments after checking for signs of variability and binarity, thus making a step further in 
mapping the top of the Cep instability strip. 

Methods. We compare the newly obtained high-resolution spectra with synthetic spectra of late 0-type and early B-type stars com- 
puted on a grid of stellar parameters. We match the results with evolutionary tracks to estimate stellar parameters. We use various time 
series analysis tools to explore the nature of the variations present in the light curve. Additional calculations are carried out based on 
distance and historical position measurements of the components to impose constraints on the binary orbit. 

Results. We find that HD 51756 is a wide binary with both a slow (vsin/ ^ 28kms~^) and a fast (vsin/ ^ 170kms~^) early-B 
rotator whose atmospheric parameters are similar (T^q ^ 30000 K and logg ^ 3.75). We are unable to detect pulsation in any of the 
components, and we interpret the harmonic structure in the frequency spectrum as sign of rotational modulation, which is compatible 
with the observed and deduced stellar parameters of both components. 

Conclusions. The non-detection of pulsation modes provides a feedback on the theoretical treatment, given that non-adiabatic compu- 
tations applied to appropriate stellar models predict the excitation of both pressure and gravity modes for the fundamental parameters 
of this star. 

Key words. Stars: variables: general - Stars: early-type - Stars: oscillations - Stars: individual: HD 51756 - Stars: rotation - Stars: 
binaries: spectroscopic 



1. Introduction 

One of the major observing programmes for the asteroseismol- 
ogy channel of the CoRoT satellite ( [Baglin et al.||2006) , is to 
derive details of the internal physics of st ars within the insta - 
bility strip of early-type B pulsators (e.g., Michel et al.||2006 



A major motivation of this programme was to understand the 
excitation problems reported for some of the observed modes 
in several prototypical bright j3 Cep stars by groups using inde- 
pendent stel lar evolution a nd pulsa tion codes (e.g .,iPamy atnykh 
leTal. 2004; Ausseloos et al. 2004; Handler et al. 2006; Desmet 



The CoRoT space mission was developed and is operated by the 



French space agency CNES, with participation of ESAs RSSD and 



Science Programmes, Austria, Belgium, Brazil, Germany, and Spain. 

Based on data gathered with Coralie installed on the 1.2 meter 
Eule r telescope at La Silla, Chile; and Harps installed on the 3.6 meter 
ESO telescope (ESO Large Programme 182.D-0356) at La Silla, Chile 
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|et al.|[2 QQ9] ), while other class members show oscillations as 
predicted by theory (|Aerts et al. 2003; Mazu mdar et al.||2QQ6| 
iBriquet et al.|2QQ7b Handl er 2009 ). The goal was also to study 
stars of spectral type O or BO, which are higher up in the pre- 
dicted 13 Cep strip than most class members observed from the 
ground prior to CoRoT, to see if the emptiness of that part of 
the strip is removed when going from mmag to yumag precision 
in photometric data. A third motivation was to obtain data of 
higher precision than from the ground and uninterrupted, in or- 
der to establish once and for all the firm detection of non-radial 
gravity modes seemingly present in ground-based spect roscopy 
of rapidly rotating Be stars (e.g., Rivinius et al.||2Q03) , and to 
interpret those in terms of oscillations of rapid rotators. A num- 
ber of early-type B stars were thus selected to fulfill these goals, 
keeping in mind the strong pointing restrictions of the satellite 
as well as the brightness limitations (e.g., |Auvergne et al.|2Q09 



The star HD 51756 (B0.5 IV) studied in this paper was selected 
as target in this context, as well as 6 O stars measured during a 
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short run. All further planned B targets have later spectral type 
due to pointing restrictions of the satellite. 

The slowly rotating B stars observed by CoRoT's astero- 
seismology CCDs so far, have led to some remarkable results. 
The only accessible known ySCep star, HD 180642 (B1.5 II-III), 
turned out to have numerous self-excited modes rather than one 
dominant radial mode as w as thought prior to CoRoT ( Degroote 
[etal. 2 009 ; Briq uet et al.|2 009) and many of the observed modes 
are resonantly coupled. Moreover, the star was found to hav e 
modes of a stochastic nature as well ( |Belkacem et al.||2009| ), 
which was never found before in a ySCep pulsator and which 
requires the presence of a surface convection zone capable to 
excite such modes. Stochastically excited modes exhibiting a 
constant frequency spacing were, surprisingly, detected in the 
08.5 V star HD 46149 (Degr oote et"aLl|201Qb| ). On the other 
hand, the 09V star HD 46202 turns out to be the most mas- 
sive star known to date with ySCep-like pulsation frequencies, 
but for which present excitation computati ons fail to predict the 
observed oscillations ( Briquet et al. |2010 ). Further, at the cool 
border of thcySCep instability strip, the B3 V star HD 50230 re- 
vealed oscillations of hybrid nature, i.e., the presence of both 
gravity modes as expected in slowly pulsating B stars and pres- 
sure modes as in the classical ySCep stars. The gravity modes 
of this star displayed a period spacing and a periodic deviation 
thereof, which allowed one to prove the existence of a chemi- 
cally inhomogeneous zone around the fully mixed core which 
is not yet tak en into account in current B star models ( [Degroote 
let al.|2010aD . 

As for the rapidly rotating early-Be stars, HD 49330 
(BO. 5 IVe) turned out to be a clear case where multiple mode 
beating was observed in real time and this phenomenon was 
found to be connected with an outburst feedi ng the circumstellar 
disk ( ,Huat eTaLl|20Q9l [Floquet et al.|[2Q09l). This su pports the 
early findings of fBaade et al.| ( |1988| ) and |Rivinius e t al. (2001) 
that the outbursts detected in the B2e star [i Cen have a pulsa- 
tional origin and are connected with multimode beating. This 
outburst phenomenon was so far not found for the other Be stars 
which were monitored i ntensively by CoR oT, but those targets 
are of later spectral ty pe ( |Neiner et al.'[2009l HD 181231, B5 IVe; 



Gutierrez-Soto et al.|2 009 , HD 175869, B8 Ille). There is yet no 



acceptable seismic modeling of the interior structure of a Be star, 
due to the mathematical complexity of the treatment of stellar 
oscillations in a deformed star. 

It is obvious from these recent findings that the diversity of 
stellar oscillations in and near the p Cep instability strip is much 
larger than anticipated before the CoRoT mission. This must im- 
ply that details in the internal physics of the various stars, such 
as their internal rotation, mixing, settling and radiative levita- 
tion due to atomic diff'usion, etc., may be diff'erent. Some prob- 
lems may also be solved by increasing the opacity in the ex- 
citation layers. In an attempt to increase the number of well 
studied early-B stars, HD 51756 (BO. 5 IV) was monitored by 
CoRoT's seismology CCDs and studied in parallel with ground- 
based spectroscopy. 

2. Fundamental parameters 

2.1. Prior to our study 

HD 51756 - a bright (V^ag = 7.23) B0.5 IV field-star in the con- 
stellation of Monoceros - was first listed among early-type high- 
luminosity objects by Miinch (1951), being classified as a B3 
star from objective-prism plates. The determination of the lu- 
minosity class and a more precise spectral type of BO. 5 IV was 



done by Morgan et al. (1955]), which is generally used since 
then. Although there are exceptions: a BO III classification by 
IDworetsky et aL]fl982 ), a B 1 V by[ Abt| ( |20Q8] ), and a B 1/2 lb by 



Wri ght et al. ( |2003| r^ there seems to be a quite good agreement 



on the spectral type, but not the luminosity class. 

Measurements by |Bumhain| ( |1875] ) showed that HD 51756 
is a triple system (catalogued as BU327, but widely known 
as WDS J06585-0301), consisting of two almost equally bright 
(Am ^ 0.3 mag) close components with a separation of p 0.7'', 
and a fainter star further away (Am ^ 4 mag, p ^ 13.7''). 
Unfortunately, only one historical radial velocity (RV) measure- 
ment (Vrad = 25 + 5 km s~^ ) can be found in th e literature ( Wilson 
1953 ). Searching for Cep stars, |Hill| ( [l967| ) took the first pho- 
tometric time-series and concluded that this field-star was con- 
stant. Among other parameters, he also measured the projected 
rotational velocity to be vsin/ = 30kms"^ There is also a 



measurement from |Glebocki & Stawikowski| pOOO| ) claiming a 
vsin / = 40kms~^ 

Further information we found concerns y^ag = 7.23, /^mag = 
2.569, Mv = -4.1, the colours of{B-V) = -0.10 and (U-B) = 
-0.93, and the colour excess of E(B - V) = 0.19. In addition 
to th ese, (U\^- V) = -2.60 and E(UV) = 1.07 was measured 



by Weber et al. (1971). A distance modulus of Vq - My = 10.7 



is given by | Vogt (1976), which leads to a distance of 1.38 kpc 



(1.04 kpc with extinction taken into account). There are simi- 
lar rn easurements of 1.3 ± 0.5 kpc ( Avedisova & Kondratenko| 
1984 for the BFS 55 H n regioij^ which is claimed to be ionised 
by the nea rby multiple subgian t, HD 51756), 1.6 kpc ([Savage 



etal. 



1985 ) and 2.156+|?-^^^ kpc (Kaltcheva & Hilditch 



2000 



with 



= -^.i)15 from uv byfi photometry). Ihe original flipparcos 
(Ferr yman et al. 1 1 997 ) p arallax is 0.56 + 4.8 mas, while the value 
from the new reduction ( [van Leeuwen|2007) is - 1 .94 + 1 . 1 3 mas. 
Unfortunately none of these lead to a useful distance estimate. 

Empirical temperature calibrations (jGulati et al.|1989|) gave 
values of log T^^ = 4.522 ± 0.051 K, log T^^ = 4.532 ± 0.042 K, 
and logTeff = 4.530 ± 0.045 K, derived from QUV, (ml965 - 
y)o, and (B - V)o, respectively (corresponding to - with errors 
taken into account - an eff'ective temperature range of 29 500 - 
37 500 K). 

2.2. New spectroscopy 

In the framework of ground based preparatory and follow-up ob- 
servations for the CoRoT space mission, high resolution and 
medium to high signal to noise spectra were ta ken with the 
CoRALiE (Baranne et al. 1996, Quel oz et al.|[20QT| R ^ 50000) 
and Harps (Mayor et al. 2003, 7? ^ 80000) instruments (78 and 
15 measurements respectively - see Table[T]for a summary of the 
spectroscopic observ ations). In addition to these, one spectrum 
gathered with Feros ( [Kaufer et al.|1999| 7? ^ 48 000) was taken 
from the GAUDI archive (Solano et al.|2005] ). 

We performed a careful normalisation of each spectrum, us- 
ing cubic splines which were fitted through some tens of points 
at fixed wavelengths, where the continuum was known to be 
free of spectral lines and telluric features. The final selection 
of these nodal points was done after several tests and quality 
checks, making sure that, e.g., the wings of the Balmer-lines are 
not aff'ected, but all artificial features (like periodic waves in the 
continuum) are corrected by the process. Cosmic removal via 



^ the UV filter is cente red at 1500A 
for the same region Fich & Blitz||l984| ) found a kinematic distance 



of 1.89 + 0.73 kpc determined by using CO velocities of the associated 
molecular clouds and the CO rotation curve of the outer galaxy 
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Table 1. Logbook of the new spectroscopic observations of HD 51756 obtained between January 2009 and December 2009 - 
grouped by observing runs. 



Instrument 


N 


HJD begin 


HJD end 


(SNR> 


SNR-range 


T 

-■-exp 


R 


CORALIE 


53 


2454842 


2454852 


124 


[97, 142] 


1500 


50000 


CORALIE 


25 


2454912 


2454922 


130 


[99, 160] 


[1500, 2000] 


50000 


Harps 


15 


2455174 


2455195 


293 


[265, 330] 


[450, 840] 


80000 



Notes. For each observing run, the instrument, the number of spectra N, the HJD, the average signal-to-noise ratio (calculated in the line free 
region of [5370A, 5400A]), the range of SNR values, the typical exposure times (in seconds), and the resolution of the spectrograph are given. 




1860 1900 1940 1980 
Date of measurement (year) 



1860 1900 1940 1980 
Date of measurement (year) 



Fig. 1. Position angle (left) and separation (right) measure- 
ments of the two bright and close components of HD51756 
(WDSJ06585-0301 AB) from The Washington Double Star 
Catalog, with a linear trend fitted to them. 



pixel-by-pixel sigma clipping and a sophisticated order-merging 
(taking into account the signal to noise values in the overlap- 
ping ranges, and correcting for the sometimes slightly diff'erent 
flux-levels of the overlapping orders) was also done by the same 
semi-automatic script. This is much faster and less subjective 
than one-by-one manual normalisation. 

Radial velocities were measured from cross-correlated pro- 
files - which were calculated from the three strongest Si iii lines 
(at 455 3 A, 4568 A, a nd 4575A) using th e least-square decon- 
volution technique of Donati et al. ( 1997| ) - by fitting a simple 
Gaussian to them. It is clear from the determined radial veloc- 
ities, that the orbital period has a timescale much longer than 
the timespan of the spectroscopic observations, and the orbit can 
not be fitted on this number of measurements. The idea of a very 
long orbital period is confirmed by position measurements cov- 
ering 122 years found in the Washington Double Star Catalog 
also (see Table ] A. l) , which show only minor change during more 
than a century in position angle (which slowly increases), and in 
separation (which slowly decreases) - see Fig.[T] 

By visual inspection of the averaged spectrum of both instru- 
ments (and the original spectra), we conclude that HD 51756 is 
a double-line spectroscopic binary (SB 2) with a slow (referred 
to as the primary from now on, as the sharp spectral features are 
connected to this star) and a fast rotator (secondary) star (see 
Fig. [2]). Seeing both components in the spectrum is not surpris- 
ing as the separation of the components is only 0.7'', which is 
smaller than the fiber aperture of Coralie (T') and Harps (V), 
and it is also below the average seeing of La Silla (0.9''). while 
the diflference in brightness of the primary and the secondary is 
small. 

Knowing the RV values corresponding to each exposure, we 
transferred all spectra to the laboratory rest frame (of the pri- 
mary). At this point, we re-run the normalisation script (this 
way the nodal points are really at the same rest wavelength), 
and calculated the average spectrum for both instruments - us- 
ing weights according to the calculated or given SNR values of 
each spectrum. The two average spectra calculated from Coralie 



and Harps are almost perfectly identical (with an SNR [5370- 
5400A] of 873 and 1016, respectively). The spectra show only 
absorption features, except the Feiii emission line at 5243A. 

Though in the case of the individual Coralie spectra, the 
SNR levels were a bit lower than one would typically have for 
such analysis, we searched for line-profile variations using the 
FAMIAS software package (jZimaj2008]), but without any posi- 
tive detection. 



2.3. Fundamental parameters 

An SB 2 binary is extremely useful when the orbit can be fitted 
and, after spectral disentangling, a full analysis can be carried 
out on both components separately. In our case, this is not possi- 
ble, due to the lack of knowledge concerning the orbital param- 
eters. 

To have a relatively fast but still accurate solution, we de- 
cided to carry out a full grid search with seven free parameters 
(Tqs, logg, and vsin/ for both components, plus the metallic- 
ity Z which was assumed to be the same for both stars), us- 
ing the OSTAR2002 jhmz & Hubeny||2QQ3] ) and BSTAR2006 
( [Lanz & Hubeny 2007| ) grids. These sets of synthetic spec- 
tra were both calculated from NLTE, plane-parallel, hydrostatic 
model atmospheres, with a microturbulent velocity of lOkms"^ 
and 2 km s"^ for the OSTAR2002 and BSTAR2006 grid, respec- 
tively. Equivalent widths of the instrumental line profiles were 
also taken into account while calculating the broadening of the 
synthetic composite spectra. Chi-square (x^) values were cal- 
culated using only regions of selected H i. He i, Si iii, and Si iv 
lines (some example fits are shown in Fig.|2]). Errors were es- 
timated by running the grid search using both the Coralie and 
Harps average spectra, using all four possible combinations of 
the two grids for the two components, and using more or fewer 
line regions for the calculation, then finally taking the most 
common parameters among the combinations with the lowest 
values as the final result, while the lowest and highest param- 
eter values close to the minima of the set the error bars. A 
summary of the results is shown in Table|2] The two stars are 
very similar, the only diff'erence being a diff'erent equatorial ve- 
locity or a diff'erent stellar inclination. The errors on the eff'ec- 
tive temperature and on the surface gravity might be slightly un- 
derestimated in case of the secondary, as this component con- 
tributes less to the overall appearance of the composite spectra 
because of the shallow lines caused by its high projected rota- 
tional velocity. The lowest x^ values were reached by using the 
OSTAR2002 grid for the primary and the BSTAR2006 grid for 
the secondary (note that, not only the microturbulent velocity, 
but also the atomic data included in these model atmosphere 
grids are slightly diff'erent). Also a slight trend in the full set 
of results suggests that the secondary might be ^ 1000 K cooler 
than the primary and it might have a slightly lower log g value 
compared to its companion star, but these deviations are within 
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Table 2. Observed and corresponding model parameters 



Rarameter 



Rrimary 



Secondary 



from spectroscopy 



logg(cgs) 
vsin/(kms~^) 
Z 



30000 + 1000 
3.75 + 0.25 
28 + 4 



30000 + 2500 
3.75 + 0.25 
170+15 



e [0.010,0.020] 



from evolutionary tracks 



M(Mq) 


18 9+9.1 


18 9+^^-^ 


R(Rq) 


9-6!^:? 


^•"-3.7 


logL(Lo) 




4 Q+0.6 


age (Myr) 




7 6+2-^ 



the given error bars. As shown below, this is in agreement with 
the observed Am = m2 - mi = M2 - Mi ^0.3 mag brightness 
difference of the components. From the definition of the mag- 
nitude scale, and the relation between eff'ective temperature and 
luminosity, we write 




4.6 4.5 4.4 4.3 

logTeff (K) 



4.6 4.5 4.4 4.3 

logTeff (K) 



Fig. 3. Overview of the main sequence evolutionary tracks of late 
0-type and early B-type stars for masses of 12 - 35 Mq (with a 
step of 0.5Ato) and a core overshoot parameter ofov = 0.2 pres- 
sure scale heights. Grid-points which fall within the error boxes 
of the primary and secondary components are shown with cyan 
and blue shading, respectively, on the log TQ^-logg (upper left), 
log Teff-log L (upper right), log Tq^-R (lower left), and log Tqq- 
age (lower right) diagrams. 



[R2 \T2 



^q(M2-Mi)/2.5 _ -^^3 



and, using this expression, we can give the upper limit of the 
diff'erence in eff'ective temperature and radii of the components. 
Assuming the same radii, we need T2 = 0.937 Ti to fulfill the 
equation, while assuming the same temperature requires R2 = 
0.877 7?i . At a Tq^ = 30 000 K this means that the secondary can 
be cooler by a maximum of ^ 2000 K, which fits the results of 
the grid search within the errors. 

After we determined the Tq^ and log g values of both compo- 
nents, we matched them w ith evolutionary track s (see Fig. |3] and 
Fig.gfor an overview, and |Briquet et al.| ( |2010| ) for a description 
of the input physics). As the lowest values indicated a mini- 
mum in between Z = 0.02 and Z = 0.01, we fixed the metallicity 
at Z = 0.014 in accordance with Przybilla et al. ( 2008) who pro- 
posed the cosmic abundance standard (CAS) from early B-type 
stars in the solar neighborhood. We choose X = 0.715 also in 
good agreement with the proposed CAS, and set the core over- 
shoot parameter at = 0.2 pressure scale heights, since aster- 
oseismic modeling results of JS Cep targets have giv en evidence 
for th e occurrence of core overshooting of that order ( Aerts et al. 
2010| ). From the model points which lie within the error box set 
by the Tq^ and log g determination, we estimated the stellar pa- 
rameters of the components. These values and their uncertainties 
are listed in TableO 



2.4. Constraints on the orbit 

Using the masses from the evolutionary tracks (see Table|2]), the 
distance- values (see Sect. |2.1| ), and t he visual separation from the 



position measurements (Table A. 1 ), we calculate the expected 



orbital period of the AB components of the system to see if 
the slow change in position angle and separation is compatible 
with it. Assuming the simplest case, a circular, pole-on orbit, the 
semi-major axis (hence the radius in this simplified case) of the 
orbit is a = 0.0056^qqq3q pc (using a dista nce o f d = I.^^q^ kpc 








4.65 4.60 4.55 4.50 4.45 4.40 4.35 4.30 4.25 

logTeff (K) 

Fig. 4. log Teff-log g diagram of late 0-type and early B-type 
stars. The spectroscopic parameters and corresponding error es- 
timations are indicated by thick and thin crosses for the primary 
and secondary component, respectively. The main sequence evo- 
lutionary tracks are for masses of 12 - 35Mq (with a step 
of lAlo) and a core overshoot parameter ofov = 0.2 pressure 
scale heights, and their shading and thickness indicates the cor- 
responding radius values. Various jSCcp instability strips for 
Z = 0.020 (continuous lines) and Z = 0.010 (dashed lines), 
both for four diff'erent combinations of metal mixtures and opac- 
ity computations (OPAL GN93, OP GN93, OPAL AGS05+Ne, 
and OP AGS05+Ne) are plotted with red, b lue, black, and green 
colours, respectively. These are taken from |Miglio et al. ( |2007 ) 
who computed them for masses up to 18 At©. 



of p = 0.68 ± 0.2'' calculated from the most recent and precise 
interferometric measurements between 1994 and 1997). Using 



^G(Ma + Mb) 



representing the values mentioned in Sect. 2.1 and a separation 



we obtain an orbital period of P = 6385^^^^^ year - with the 
errors covering all uncertainties. These values suggest a change 
of 6.9^4^j-^ degrees in position angle during 122 years, which fits 
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1.0 



S 0.9 
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^ 1-00 
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^ 0.85 
0.80 
0.75 



4840 
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1 


1 1 1 1 1 



1.10 



4003 4006 4009 
A (A) 
Si IV 



4012 4015 



1.05 - 



% 1.00 



I 0.95 



4550 4555 4560 4565 4570 4575 
A (A) 



0.90 - 
0.85 



1 r 




J L 



4208 4210 4212 
A (A) 



4214 4216 



Fig. 2. Fits to the averaged Harps spectrum (black). The parameters are listed in Tabled The fit to the primary companion is shown 
in solid grey, the fit to the secondary in dashed grey, and the combined fit in red. While these synthetic spectra were calculated for 
Z = 0.010, we also plot the combined spectrum corresponding to Z = 0.020 with a solid blue line to show the eff'ect of the change in 
metallicity. The eff'ect of the fast rotator companion is clearly visible in the He i line at 47 13 A {lower left panel). We had difficulties 
in fitting the metal lines equally well with any given setting. Some regions suggest higher (e.g., the Sim lines at 4553A, 4568A, and 
4575A {lower middle)), while some others suggest lower metallicity values (e.g., the Siiv line at 4212A {lower right)). 



the observations (showing a change of ^ 6.3°). Assuming an 
edge-on orbit, we derive an upper limit of the semi-amplitude of 
the radial velocity; Arv = 5.4^^ g kms~^ This - combined with 
the shallow profiles of the secondary due to its high projected 
rotational velocity and the relatively short timespan of the mea- 
surements (compared to the orbital period) - explains why there 
was no (measurable) wavelength- shift between the lines of the 
two components in the spectra. 



3. The CoRoT light curve 

3.1. Frequency analysis 

HD 51756 was observed by the CoRoT satellite during the sec- 
ond long run of the mission (LRa02) as part of the asteroseismol- 
ogy programme, while monitoring a field in the direction of the 
Galactic anticenter from 2454784.491822 HJD for almost 115 
days (from 13 November 2008 to 8 March 2009). All flagged 
observations were removed from the light curve, leaving us with 
278 131 measurements (resulting in a duty cycle of ^ 90%). The 
decreasing trend of the CoRoT light curve has a well known in- 
strumental origin ( Auvergne et al. |2009 ), hence it was also re- 
moved by dividing by a linear fit. As after this step there were 
no clear jumps or trends visible anymore, we used the resulting 
dataset in our analysis. 



The first visual inspection of the light curve showed a well 
traceable repetitive pattern over almost all the timespan of the 
observation, with a length of ^ 2 days (see Fig. [5]). The shape 
and features of the light curve show a noteworthy resemblance 
to the 08.5 V star HD 46149 observed during the SRa02 shor t 
run of CoRoT and analysed in detail by Degroote et al.|( 2QTQb| ). 



We performed an iterative prewhitening procedure (see e.g. 
Degroote et al.||2009|) us ing the traditional linear Scargle peri- 
odogram ( |Scargle|| 1 982| ) . This resulted in a list of amplitudes 
{Aj), frequencies {fj), and phases {6j), by which the light curve 
can be modelled via nj frequencies in the well-known form of 



rif 

F{ti) = c + Yj^j sin[27r(//,- + Oj)]. 

We consider a peak sig nificant if it exceeds an amplitude signal- 
to-noise ratio of 4 (see Breger et al.|1993 ). The noise level was 
calculated as the average amplitude - before prewhitening - 
in a 3 d"^ interval centered on the frequency of interest. Only 
27 peaks of the Scargle periodogram (Fig.|5]) met the signal-to- 
noise criterion, and using only the significant frequencies listed 
in Table ] A. 2| we are unable to model the light curve on a satis- 
factory level. Even using all the 93 peaks which have an SNR 
above 3, the constructed model lacks the characteristics of the 
original dataset, leaving most of the non- sinusoidal sharp fea- 
tures in the residual light curve. The noise level in the original 
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Fig. 5. (upper panel) Part of the reduced CoRoT light curve, showing distinctive sharp recurring features every 1.9 days. This can 
be easily seen around 3255 CoRoT JD, where the occurrence of features causing the highest amplitude variations is clearly periodic. 
(middle panel) The Scargle periodogram of the full CoRoT light curve showing a series of harmonics at integer multiple frequencies 
of the main peak at /o = 0.52687 ± 0.00004 d"^ (lower panel) 0p/)M-statistics of the light curve, showing the same structure as the 
Scargle periodogram. 



periodogram around 1.5-5-10 d"^ is 20.7-7.8-3.7 ppm, respec- 
tively, while it is 15.7-7.1-3.5 ppm in the periodogram of the 
residual light curve (which was created by subtracting a model 
with the parameters listed in Table |A.2| ). 

For comparison - and because of the non- sinusoidal shape 
of the signal - we also used the phase dispersion minimisation 
(RDM) method of |Stellingwerf| ( |1978| ) to calculate the &pdm- 
statistics, which shows the same global structure and features as 
the Scargle periodogram (see Fig. [5]). 

To check the stability of features in the frequency spectrum, 
time-resolved Fourier transformations were also calculated with 
different window widths. It is clearly visible from Fig. [6] that the 
amplitudes connected to the significant frequencies vary signifi- 
cantly over the span of the time series. 

Of th e 27 frequencies with SNR > 4, 11 are connected with 
/o (Table [Al2| ). All these frequency properties are a- typical for 
stellar oscillations in hot B stars and point towards another origin 
of the variability. 

3.2. Repetitive patterns 

The classical Fourier-techniques revealed a series of harmon- 
ics at integer multiple frequencies of the main peak at /o = 
0.52687 ± 0.00004 d-^ All harmonics (using a notation of 
for the ^th harmonic of /o) can be traced in the power spec- 
trum down to /q^ (with SNR > 3), plus /J^ and f^^ are also 
visible. There is a small but systematic diff'erence between the 
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Fig. 6. Short-time Fourier transformation (window = 5/o) of the 
[0.2, 6.0] d"^ region of the light curve, showing that the am- 
plitudes connected to diff'erent harmonics of the main peak at 
/o = 0.52687 ± 0.00004 d"^ (indicated as red dashed lines) vary 
significantly over time. 



observed and the exact nfo (n e Z) values. Fixing the fre- 
quency values of the harmonics to nfo and carrying out a non- 
linear least- squares fitting procedure with these peaks gives a 
worse model, with remaining higher residual values caused by 
un-prewhitened peaks right next to the removed fixed frequency 
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Fig. 7. Autocorrelation function of the CoRoT light curve (upper 
left), of the power spectrum (upper right), of the light curve cut 
into two and three pieces (lower left), and of the power spectra 
which were calculated from the pieces (two and three, plotted 
with thicker and thinner lines, respectively) of the light curve 
(lower right). 



trum is a sign of rotational modulation. To carry out this test, 
we calculated the inclination values which correspond to a ro- 
tational period of l//o (1.898 d) within the error bars of the R 
and V sin / values (from Table[2]), and we checked if the equato- 
rial rotational speed at this set-up was below the critical veloc- 
ity for the given R, M combination. If we suppose that the pri- 
mary is responsible for the repetitive features in the light curve, 
then its inclination has to be /a = 6.3^3 ^ degrees (Fig.pb. On 
the other hand, assuming that the variations are connected to 
the secondary, then we conclude an /b = 41.6^21? ^^S^^^^ 
we end up with an equator-on solution already slightly before 
reaching the lowest possible R value. Taking into account only 
the most probable R values we end up with /a = ^-^-o 9 degrees 
and /b = 41.6^4 4 degrees. 

As it is clearly visible on Fig. [8] these setups may result in 
exceeding the break-up velocity if the radius of the given com- 
ponent is in the very upper part of the error box. Lowering the 
mass, the radius limit also drops. We can conclude that most of 
the possible configurations are indeed compatible with our hy- 
pothesis of a rotational origin of the variability. 



harmonics. The observed frequency values are always higher 
than the ones used in the fitting (fixed at integer multiples of 
/o), and this diff'erence increases with frequency (or n), starting 
with A/q = 2/0 - /o ~ -0.008 d"^ and continuously growing to 
A/20 = 20 fo - f^^ ~ -0.05 d-\ This might be a sign of rota- 
tional modulation with the rotational period of one of the com- 
ponents. If diff'erential rotation were present, it could complicate 
the structure of the peaks and it might explain in part the diff'er- 
ences between the observed peak values and exact multiples of 
fo. 

Another approach to search for spacings and recurrent fea- 
tures both in the light curve and the power spectrum is to use 
the autocorrelation function. We calculated the autocorrelation 
of the light curve, and the power spectrum (here the square of 
the Scargle periodogram). These tests were repeated after the 
dataset was cut into two, and three pieces, respectively. The re- 
sults displayed on Fig.|7]show clear signs corresponding to the 
repetitive nature of the variations in the light curve. This signa- 
ture stays quite clear even for shorter subsets of the data, but 
gets less significant in case of the power spectrum calculated 
from these data sets. Furthermore, the location of the first peak 
of the autocorrelation function stays almost perfectly the same 
even for three subsets, but this stability can not be seen in case 
of the autocorrelation of the power spectra of the subsets. 

3.3. The origin of repetitive variations 

A possible source of the variability might be atmospheric fea- 
tures, like spots or chemical inhomogeneities on the surface. In 
that case the observed brightness variations can be explained by 
the appearance and disappearance of such regions, either be- 
cause of repeated creation or destruction or because of the ro- 
tation of the star. Here, we are in favour of the second option, 
because of the signal's self- similarity, and because - as we show 
it below - it is consistent with the determined v sin / value. For 
further details on how these features can actually produce the 
observed cusp-like peaks, we refer to |DegrooTe| ( |20 1 0| ) . 

Knowing the projected rotational velocities of both compo- 
nents, we test if it is physically possible that one of the com- 
ponents has a rotational period which is compatible with our 
hypothesis that the series of harmonics in the frequency spec- 



3.4. Inferences from the non-detection of pulsation 

For all the models passing through the error boxes of the primary 
and secondary indicated in Fig.|5j we considered the frequencies 
of the fundamental radial mode and of the few lowest radial over- 
tones, as well as the lowest-order p and g axisymmetric modes 
of degree 1 and 2. None of those eigenfrequencies come close 
to fo. It is thus very unlikely that fo, with so many harmonics 
present in the light curve, results from a pulsation mode, as one 
does not expect a high-overtone mode or a high-degree mode to 
show so many harmonics of its frequency (e.g., D egroote et al. 
|2Q09i ). 

We made non-adiabatic comp utations (with the code MAD 
developed by [Dupret et al.||2002| ) for the mentioned models to 
check which modes are predicted to be excited by the present 
input physics and theory of mode excitation (see Fig. [9]). We need 
to understand why the star does not oscillate, as the theoretical 
calculations predict p modes to be excited throughout the error 
box of HD 51756, except the very bottom left region. Typical 
frequencies are 3.5 - 6.5 d"^ around the center of our error box. 
This range extends until 1.9 d"^ and 9.4 d"^ towards the lower 
and higher gravity values, respectively. Moreover, g modes are 
also expected to be excited in some of the higher mass models, 
with frequencies near 0.1 d"^ 

Finally, we should remember that CoRoT measured the inte- 
grated flux from two stars, therefore the observed amplitude of 
the rotational modulation is damped, and this would be true also 
for the amplitudes of the pulsational modes, if any. 

4. Conclusions 

We analysed high-resolution spectroscopy and high-quality 
space based photometry of HD 51756. The former - via syn- 
thetic spectrum fitting and matching the results with evolution- 
ary tracks - enabled us to derive the fundamental stellar param- 
eters of both components of the underlying double lined spec- 
troscopic binary, while the latter led us to conclude that none of 
these stars show pulsations. We interpret the repetitive patterns 
of the light curve, the observed harmonic structure in the fre- 
quency spectrum, and the strong, stable peaks in the autocorrela- 
tion functions as signs for rotational modulation. This hypothesis 
is compatible with the observed and deduced stellar parameters. 
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Fig. 8. Possible Vrot combinations in agreement with the observed vsin/ = 28kms"^ of the primary (left) and the vsin/ = 
170 km s"^ of the secondary component (right). Plots show (via colours and thin black curved contour lines) rotation periods 
calculated from the observed projected rotational velocities at given /, R values. The thick contour line corresponds to the main peak 
at 1.898 d from the frequency spectrum, and the dashed, continuous, and dot-dashed red lines show the break-up velocity (with 
speeds rising towards the left) for Almax, for At, and for Atmin, respectively. The thick contour line corresponds with the parameter- 
combinations which fit the hypothesis in Sect. |3.3| and the white dot shows the location of the most probable R (and thus /) value 
from Tabled 
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Fig. 9. log Teff-log g diagrams showing an overview of models where diff'erent modes are predicted to be excited in our non-adiabatic 
computations for the error box of the components of HD 51756. Each filled circle represents a model along the evolutionary tracks 
where modes are predicted to be excited. The size of the symbols refers to the number of excited modes (1-4 per mode type), while 
the colour corresponds to the frequency diff'erence between /o and the closest mode in the model. The spectroscopic parameters 
and corresponding error estimations of HD 51756 are indicated by thick and thin crosses for the primary and secondary component, 
respectively. 
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The absence of oscillations in the components of HD 51756 
requires an explanation. It was recently found by Balona et al. (in 
preparation) that some stars observed with the Kepler space mis- 
sion, which are also monitored at yumag precision and situated in 
the p Cep instability strip, do not pulsate either. A mixture of 
non-pulsators and pulsators with or without rotational modula- 
tion was already established from ground-based data for stars 
in the g-mode instability strip of the cooler slowly pulsating B 
stars (e.g., |Briquet et al.||2004| ). It was interpreted as due to a 
difference in magnetic field strength and age, the Bp stars being 
younger and having stronger fields than the pulsators, although 
this int erpretation is still a matter of debate (e.g.. Briquet et aLj 
2Q07a Silvester et al. 200 9 ). As far as we are aware, no magnetic 
field measurements are available for HD 51756. 
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Appendix A: Tables 

Table ] A. 1 1 contains all the position measurements of HD51756 
(catalogued as BU327, but widely known as WDSJ06585- 
0301) found in the Washington Double Star Catalog. Error esti- 
mates - when available - can be found in the original references. 

Table ] A.2| lists the Fourier parameters of significant peaks in 
the periodogram. 



Table A.l. Position measurements of HD 51756 AB 



Date 


on 


pD 


N 


Method"" 


Reference^ 


1875.06 


90 


1 


1 


A 


1 


1876.83 


100.8 


0.96 


2 


A 


D1883 


1878.12 


97.2 


n/a 


1 


A 


2 


1881.05 


96.1 


0.75 


1 


A 


Bul883 


1888.10 


95.8 


0.73 


2 


A 


3 


1888.222 


95.5 


0.84 


4 


B 


4 


1892.08 


95.8 


0.79 


3 


A 


Bul894 


1894.40 


98.1 


0.64 


5 


A 


Spl909 


1899.01 


93.6 


0.59 


1 


A 


Bul900 


1902.544 


93.9 


0.68 


6 


A 


5 


1903.14 


100.0 


0.64 


2 


A 


9 


1906.16 


95.9 


0.87 


3 


A 


6 


1906.17 


101.5 


0.81 


1 


A 


Frml907 


1909.06 


100.4 


0.90 


1 


A 


7 


1909.06 


101.3 


0.88 


1 


A 


7 


1909.13 


96.9 


0.70 


2 


A 


Wzl912 


1911.15 


91.7 


0.56 


2 


A 


8 


1920.09 


98.1 


0.77 


1 


A 


10 


1924.35 


99.4 


0.83 


3 


A 


10 


1924.79 


97.5 


0.94 


1 


A 


10 


1930.57 


98.3 


0.71 


2 


A 


Bonl938 


1930.84 


96.2 


0.80 


3 


A 


11 


1935.43 


102.0 


0.78 


4 


A 


12 


1942.88 


95.1 


0.62 


3 


A 


13 


1945.31 


103.6 


0.68 


2 


B 


VBsl954 


1947.17 


102.8 


0.76 


2 


A 


VBsl954 


1953.22 


98.5 


0.81 


6 


A 


Rabl961b 


1956.11 


98.2 


0.74 


2 


A 


14 


1961.04 


101.6 


0.67 


4 


A 


15 


1965.034 


100.1 


0.69 


4 


A 


16 


1966.182 


101.4 


0.58 


3 


A 


17 


1968.045 


101.4 


0.66 


3 


A 


18 


1974.717 


101.1 


0.59 


3 


A 


19c 


1975.522 


101.6 


0.55 


3 


B 


Wakl985 


1991.25 


102. 


0.71 




T 


20 


1991.92 


104.0 


0.68 


I 


T 


21 


1994.129 


101.2 


0.67 


I 


s 


22 


1994.194 


99.7 


0.63 


I 


s 


22 


1994.194 


100.5 


0.68 


1 


S 


22 


1995.188 


100.6 


0.66 




S 


22 


1995.199 


103.0 


0.68 




S 


22 


1995.202 


102.1 


0.70 




S 


22 


1995.202 


101.8 


0.67 




S 


22 


1995.205 


101.5 


0.69 




S 


22 


1995.207 


101.9 


0.63 




S 


22 


1996.8960 


101.5 


0.695 




S 


23 


1997.184 


101.5 


0.67 




S 


24 


1997.184 


102.2 


0.69 




S 


24 



Notes. Methods of observations. A: Refractor, micrometer; 
B: Reflector, micrometer; T: Hipparcos or Tycho observation; S: 
Speckle interferometry.*^^^ For special references - which are not 
in the SAO/NASA Astrophysics Data System - see the References 
and discoverer codes of The Washington Double Star Catalog at 
|http://ad.usno.navy.mil/wds/Webtextfiles/wdsnewref.txt'^^ These obser- 
vations were made after the closing of the referred article, but before the 
retirement of Alan Behall. Observing t echniques, ins trumentation and 
particulars are the same as described by [Behall 1(1976] ) 
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Table A.2. Table of Fourier parameters (frequencies (fj), amplitudes (Aj), and phases (6j)) of peaks with a signal-to-noise ratio 
(SNR) > 4, plus clear harmonics (f^) of the main peak at /o = 0.52687 ± 0.00004 d"^ down to SNR > 3. SNR is calculated over 
3 d"^ in the periodogram. 



/(d"0 


6/(d 


A (ppm) 


6a (ppm) 


6 (2;7r/rad) 


ee (2;r/rad) 


SNR 


note 


0.009071 


0.000035 


167.4 


1.1 


0.2563 


0.0073 


10.4 




0.017911 


0.000030 


183.8 


1.1 


-0.1304 


0.0063 


12.2 




0.033302 


0.000042 


137.5 


1.1 


-0.0228 


0.0087 


8.5 




0.040540 


0.000060 


80.1 


1.0 


0.3908 


0.0124 


6.0 




0.062041 


0.000057 


89.8 


1.0 


0.2046 


0.0118 


6.4 




0.098717 


0.000058 


102.2 


1.0 


-0.0143 


0.0121 


6.2 




0.119596 


0.000080 


58.2 


0.9 


-0.0415 


0.0165 


4.5 




0.137555 


0.000083 


59.1 


0.9 


0.4332 


0.0173 


4.3 




0.253269 


r\f\f\r\'~i c 

0.000075 


60.2 


0.9 


f\ ^ A 0'~l 

0.3487 


0.0156 


4.7 




U.JUOOZD 


U.UUUU/1 


dJ.9 


1 rv 
l.U 


A '2 1 OA 


A Al /I 

U.U140 


C 1 

J.l 




\jA 1 lylK) 


r\ r\f\f\r\n i 
U.UUUU/ i 


DJ.J 


l.U 


-U.Z/J / 


A Al An 

U.Ui4/ 


< A 




U.4oJ JO / 


U.UUUUo i 


Oy.J 


n Q 




A Al (iO 

U.UiDo 


A A 
4.4 




n ZLQ^^QZLQ 


U.wUUUoo 


^ / .o 


n Q 


n 1 Q44 
w. 1 ^^-H- 


n ni 89 

w.wi oZ 


4 1 
4. 1 




U. JZDoDJ 


w.wwwwH-z 


1 94 1 


1 
1 .w 


^8Q7 


0088 
w.wwoo 


8 S 

O.J 


JO 




00007 S 

W.WVJWW / J 


U 1 . J 


Q 


-0 078^ 

-w.W / OD 


01 

w.WU J 


4 8 




551 246 


000060 

W. WW WW WW 


69.7 


1.0 


0.0598 


0.0124 


6.0 




1 061 ^Q6 


000058 

W. WWWWJ o 


94.0 


1.0 


-0.3826 


0.0120 


6.2 


h 




000050 

W. WWWWJ w 


101.1 


1.0 


0.3509 


0.0104 


7.2 




2.122824 


00008 1 

W.WWWWO X 


58.3 


0.9 


-0.3705 


0.0167 


4.4 


h 


2 65SQ12 


000106 

w. www ± WW 


39.1 


0.9 


-0.0592 


0.0219 


3.4 




s 1 81 son 


0000Q8 

w.wwwwvo 


42.5 


0.9 


0.2905 


0.0203 


3.8 


? 


71 8828 


0000Q4 

W. WWWW2/'-T 


45.3 


0.9 


-0.1821 


0.0195 


4.3 


H-.ZH-J iUi 


U.wUUi iU 


^1 9 


U.o 


n 1 ^1 Q 


n 0990 
w.wZZ>' 


4 1 
4. 1 


Jo 


4.769149 


0.000094 


45.5 


0.9 


0.1209 


0.0195 


5.5 


fo 


5.297008 


0.000119 


31.9 


0.8 


-0.1407 


0.0248 


4.6 


flO 
JQ 


5.809200 


0.000188 


17.3 


0.6 


0.0690 


0.0390 


3.1 


Jo 


6.354499 


0.000185 


16.4 


0.6 


-0.3715 


0.0384 


3.2 


Jo 


6.894291 


0.000213 


13.5 


0.6 


0.4758 


0.0442 


3.1 


fU 

Jo 


7.421695 


0.000163 


19.8 


0.7 


-0.2722 


0.0338 


4.2 


fl4 
^0 


7.949167 


0.000198 


14.9 


0.6 


-0.1871 


0.0411 


3.6 


Jo 


9.009843 


0.000178 


17.3 


0.6 


-0.4902 


0.0369 


4.6 


r\l 

Jo 


10.583197 


0.000140 


24.9 


0.7 


-0.3271 


0.0291 


7.4 


f20 

Jo 


13.969336 


0.000292 


9.0 


0.5 


0.1131 


0.0606 


4.0 


a 



Notes. The group of low frequency values show up because of noise and not completely removed trends. Several frequencies close to /o can be 
explained by the fact that neither the frequency nor the amplitude of /o is stable, as it can be clearly seen on Fig.[6] Peak corresponding to the 
orbital period of the CoRoT satellite. 



